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ABSTRACT 



Aims. We investigate the impact of realistic three-dimensional (3D) hydrodynamical model atmospheres of red giant stars at different 
metallicities on the formation of spectral lines of a number of ions and molecules. 
Or Methods. We carry out realistic, ab initio, 3D, hydrodynamical simulations of surface convection at the surface of red giant stars with varying 
q ■ effective temperatures and metallicities. We use the convection simulations as time-dependent hydrodynamical model stellar atmospheres to 
5_l ■ calculate spectral lines for a number of atomic (Lii, Oi, Nai, Mgi, Cai, Fei, and Fen) and molecular (CH, NH, and OH) lines under the 
(y-j assumption of local thermodynamic equilibrium (LTE). We carry out a differential comparison of the line strengths computed in 3D with the 
CS • results of analogous line formation calculations for classical, ID, hydrostatic, plane-parallel marcs model atmospheres in order to estimate the 
impact of 3D models on the derivation of elemental abundances. 
. ,— | , Results. The temperature and density inhomogeneities and correlated velocity fields in 3D models, as well as the differences between the mean 
■ 3D stratifications and corresponding ID model atmospheres significantly affect the predicted strengths of spectral lines. Under the assumption 
5— i of LTE, the low atmospheric temperatures encountered in 3D model atmospheres of very metal-poor giant stars cause spectral lines from 
5^ i neutral species and molecules to appear stronger than within the framework of ID models. As a consequence, elemental abundances derived 
from these lines using 3D models are significantly lower than according to ID analyses. In particular, the differences between 3D and ID 
abundances of C, N, and O derived from CH, NH, and OH weak low-excitation lines are found to be in the range -0.5 dex to —1.0 dex for the 
the red giant stars at [Fe/H] = -3 considered here. At this metallicity, large negative corrections (about -0.8 dex) are also found, in LTE, for 
weak low-excitation Fei lines. We caution, however, that the neglected departures from LTE might be significant for these and other elements 
and comparable to the effects due to stellar granulation. 
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1. Introduction 

Stars in the red giant phase of their evolution are character- 
ized by large stellar diameters and increased luminosities com- 
pared with unevolved objects. Because of their high intrinsic 
brightnesses, red giants represent natural targets for a wide 
range of observational studies. Stars of this luminosity class 
are especially suitable for investigations of distant stellar sys- 
tems in the Milky Way as well as in other galaxies of the Local 
Group. Giant stars are extensively used in spectroscopic anal- 
yses for tracing elemental abundances in distant stellar pop- 
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ulations. A significant fraction of the halo and disk stars in- 
cluded in various large-scale stellar abundance analyses (e.g. 
lMcWaUametal]ll995t iRyan et al.l[l996t iFulbrightlEoool) are 
inde ed giants. More recently, the ESQ "First Stars" progra mme 
(e.g. Iffill et alJflooH ICavrel et alj|2004l ISpite et alJl2005T) has 
provided a systematic and homogeneous study of the chemical 
compositions of a large sample of extremely metal-poor giant 
as well as dwarf stars ([ Fe/H] & -2.5) from the HK survey by 
iBeers et aL dl992 , 1999|). The recently discov ered extreme halo 
star HE 0107 - 5240 dChristlieb et al.ll2002l) . one of the most 
iron-poor star known to date, is also a giant. Red giants are 
also used to study Galactic metallicity gradients in the Galactic 
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disk and open clusters dYong et al.l 12005c ICarnev et alj 12005) 
and star-to-star elemental a bundance variations within globular 
clusters (for a review, see iGratton et al. 2004 ). The results of 
these and similar observational studies provide valuable infor- 
mation about fundamental astrophysical processes and are es- 
sential for our understanding of stellar nucleosynthesis and in- 
ternal mixing as well as Galactic chemical evolution. Accurate 
determinations of elemental abundances in metal-poor stars 
and red giants in particular are crucial for deducing the phys- 
ical and evolutionary properties of the very first generation 



Table 1. Details of the 3D hydrodynamical simulations. 



of stars dWeiss et al.ll2004t llwamoto et al.ll2005t iMevnet et al 



2006) and distinguishing among the various proposed scenar 
ios of chemical enrichment of the Galaxy. 

As for other late-type stars, ordinary classical abundance 
analyses of red giants involve theoretical ID model atmo- 
spheres constructed under the assumptions of plane-parallel 
geometry (or spherical symmetry), hydrostatic equilib rium, 
and flux constancy (e.g. iGustafsson & Jorgensenlll994 ). The 
simplifying approximation of LTE is also normally adopted. 
Furthermore, ID modelling of stellar atmospheres generally re- 
lies on a rudimentary treatment of convective energy transport, 
such as the mixing length theory (Bohm-Vitense 1958) or one 
of its derivatives (e.g. ICanuto & Mazzitellil 1 1 99 lb . all depen- 
dent on a number of tunable but not necessarily physically well 
motivated free parameters. In late-type stars, however, the con- 
vective zone reaches and significantly affects the regions from 
which the stellar flux is emitted. High spatial resolution imag 



ing of the solar p hotosphere (e.g. lTitle et al.ll 1 990t ISpruit et al 



1990; Carlss on et al . 2004) reveals that the surface of the Sun is 
dominated by a distinct granulation pattern reflecting the bulk 
flows in the convective zone deeper inside. Observational di- 
agnostics of stellar granulation in general are more limited, as 
the surfaces of most other stars cannot be directly resolved. 
Nonetheless, other distinguishing signatures (e.g. wavelength 
shifts and asymmetries of line profiles) of the presence of pho- 
tospheric velocity fields and correlated temperature inhomo- 
geneities can be identified in h ig h-resolution spectra of late - 
type stars JPravinsIl 98211 1987albFlAUende Prieto et al.l2002ah . 
These observable manifestations of stellar surface convection 
immediately reveal that the assumptions of hydrostatic equi- 
librium and flux constancy in ID model atmospheres of late- 
type stars are questionable and might be inadequate for high- 
precision stellar abundance analyses. Furthermore, ID mod- 
els can neither predict the strengths and shapes of lines with- 
out invoking the ad-hoc fudge factors micro- and macro- 
turbulence to account for non-thermal Doppler broadening of 
spectral lines associated with bulk flows in the stellar atmo- 
sphere. During the past three decades, realistic 3D hydro- 
dynamical simulations of stellar surface convection have be- 
come feasible thanks to the advances in computer technology 
and the development of efficient numerical algorithms (e.g. 
Nordlundlll982t iNordlund & Dravins|[l990l; IStein & Nordlundl 
1 998i lAsplund et al.l|l999|; iFrevtag et al 1 120021; iLudwig et al.l 
2002t fcarlsson et al.l 120041 IVoglerfeooH) . At present, 3D hy- 
drodynamical simulations of surface convection have primarily 
been carried out for main sequence stars and subgiants of spec- 
tral types A to M at different metaiiicities; due to the disparate 
time scales for radiation transport and convection, convection 
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" Temporal average and standard deviation of the emergent effective 
temperatures. 

h Time span of the parts of simulations used for spectral line forma- 
tion purposes. 

simulations of red giant stars appear in general more compu- 
tationally demand i ng, but they are curre ntly being developed 
dCollet et al]l2006r.lKucinskas et alJl2006h . 

Three-dimensional simulations of stellar surface convec- 
tion can be used as time-series of hydrodynamical model 
atmospheres to study the effects of photospheric inhomo- 
geneities and velocity fields on the formation of spectral lines. 
Recent analyses based on state-of-the-art 3D time-dependent 
simulations of surface convection in the Sun, dwarfs and 
subgi ants (e.g. Asplund et al. 119991: lAsplund & Garcia Perezl 



200 1[ lAsplundll2005l and references therein) indicate that the 
structural differences between 3D hydrodynamical and ID hy- 
drostatic model stellar atmospheres can have a significant im- 
pact on the predicted strengths of synthetic spectral lines and 
hence lead to severe systematic effects on the derivation of 
elemental abundances, in particular at low metallicity. In the 
present paper, following the spirit of these previous works, we 
carry out the first 3D hydrodynamical simulations of surface 
convection in red giant stars with metaiiicities ranging from 
solar down to [Fe/H] = -3. We examine the main differences 
between the structures of the 3D hydrodynamical simulations 
and ID hydrostatic model atmospheres computed for the same 
stellar parameters and study their effect on spectral line forma- 
tion and abundance analyses of giant stars. 

2. Methods 

2.1. Convection simulations 

We use the 3D, time-depende nt, compressible, exp licit, 
radiative-hydrodynamical code bv lStein & Nordlundl d 19981) to 
construct sequences of surface convection simulations of red 
giant stars with varying effective temperatures and metaiiici- 
ties. The hydrodynamical equations of mass, momentum, and 
energy conservation are solved together with the 3D radiative 
transfer equation on a Eulerian mesh with 100x100x125 grid 
points. The physical domains of the simulations are set large 
enough to cover about ten granules simultaneously in the hori- 
zontal plane and eleven pressure scales in the vertical direction. 
In terms of continuum optical depth at A — 5000 A, the simu- 
lations extend from logrsooo $ -5 down to logT 5 ooo ~ 7. The 
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depth scales are optimized to provide the highest spatial reso- 
lution in those layers near the optical surface where the vertical 
temperature gradients are steepest. For the simulations, we em- 
ploy open boundaries vertically and periodic boundaries hor- 
izontally. At each time step the radiative transfer equation is 
solved along one vertical and eight inclined rays; th e opacities 
are grouped into four opacity bins dNordlundll 19821) and local 
thermodynamic equilibrium (LTE) without continuous scatter- 
ing terms in the source function (S v = B v ) is assumed through- 
out the calculations . The adopted equation of state comes from 
Mihal as et a i1 dl988h and accounts for the effects of ionization, 
excitation, and dissociation of 15 of the most abundant ele- 
ments, as well as the H2 and HJ molecules. Continuous opaci- 
ties c ome from the Uppsala opacity package ( Gustafsson et al.1 
19751 and subseq uent updates) and line opacity data from 



Kuruczl(ll992lll993l) . 

For the present work, we have generated two series of 
3D hydrodynamical surface convection simulations of red gi- 
ant stars with surface gravity logg = 2.2 (cgs). The first se- 
ries comprises simulations of red giants with r e ff^4750 K and 
metallicities [Fe/H] = 0, -1, -2, and -3 while the second se- 
ries corresponds to red giants in the same metallicity range but 
with somewhat higher effective temperatures (7^^5050 K). 
For all hydrodynami cal simulations, we adopt th e solar chemi- 
cal composition from lGrevesse & Sauva with the abun- 
dances of all metals scaled proportionally from solar to the rele- 
vant [Fe/H] and without ^-enhancements. The initial snapshots 
for the present red giant simulations are obtained from simula- 
tions performed at a lower numerical resolution (50x50x125), 
which have been running for sufficiently long times to cover 
several convective turn-over time-scales and allow thermal re- 
laxation to occur. The initial snapshots of the lower resolu- 
tion simulations in turn have been constructed by appropri- 
ately scaling a set of previous simulations of turn-off stars by 
Asplund & Garcia Perezl ( 2001 ) to the new stellar parameters, 
using the experience from classical ID, hydrostatic stellar mod- 
els. 

Some relevant parameters and characteristic physical quan- 
tities of the red giant convection simulations are given in 
Table Q] While the geometrical depths of the simulation do- 
mains resulting from the scaling process vary depending on 
the metallicity, all convection simulations span approximately 
five pressure scale heights below and six pressure scale heights 
above the optical surface. We deem this sufficient to realis- 
tically simulate the velocity fields and spatial structures at 
the stellar surface. It is worthwhile mentioning that, in the 
convection simulations, the entropy of the inflowing gas at 
the lower boundary replaces the effective temperature as a 



const ant and independent input parameter (IStein & Nordlund 



119981) . Consequently, the emergent effective temperatures of 
the simulations vary slightly with time, the total outgoing ra- 
diative flux being susceptible to the evolution of the surface 
granulation pattern. The actual value of the temporally aver- 
aged T e g ultimately depends on the entropy of the inflowing 
gas. Constructing a simulation with a specific value of the tem- 
poral averaged effective temperature requires a careful fine- 
tuning of the entropy parameter. As this procedure is very time- 
consuming, and as we are primarily interested in a differential 



comparison between 3D and ID models, we consider it satis- 
factory for the scope of the present paper to settle for values 
reasonably close to the targeted effective temperatures for the 
two suites of convection simulations. 

From a qualitative point of view, the atmospheric structures 
and gas flows resulting from the convection s imulations are 
fairly simila r to the ones previously found by Asp lund et alJ 
dl999h and lAsplund & Garcia Perezl d200lh for dwarfs and 
turn-off stars. The warm isentropic gas ascending from the stel- 
lar interior rapidly cools as it approaches the optical surface, 
where it also loses entropy; the cooled gas eventually turns 
over and falls back toward the interior due to negative buoy- 
ancy. The morphology and the evolution of the resultant gran- 
ulation patterns, with warm, large upflows amidst cool, narrow 
downdrafts, are essentially the same as for solar-type stars. 

We have also computed classical I D, LTE, plane-pa rallel, 



hydrostatic marcs m odel atmospheres dGustafsson et al.ll 1975 



Asplund et al. 1 19971) with identical stellar parameters, input 



data, and chemical compositions as the 3D simulations to al- 
low for a differential comparison of the two types of models 
in terms of spectral line formation and abundance analysis. For 
the marcs models presented here, we ha ve adopted the mixing - 
length theory (MLT) formulation from iHenyey et alJ d 19651) . 
with, in particular, the cumlt parameter set to 1.5, the structure 
parameter y to 3/(47r 2 ) a 0.076, and the v parameter to 0.8. We 
have not considered the effects of turbulent pressure when con- 
structing the ID models. According to our tests, the inclusion 
of a turbulent pressure term P t urb = fipu 1 , with u the turbulent 
speed and/3 = 1.5, affects the gas pressure-temperature relation 
but mainly in those layers below optical depth log T5000 - 0.4. 
In addition, calculations of synthetic Fei and OH line profiles 
including and excluding the turbulent pressure term in the ID 
models indicate that the resulting differences in Fe and O abun- 
dances determinations are typically 0.01 dex or less. 

Contrary to the convection simulations, in the ID models, 
scattering is correctly treated as such and not as true absorbtion. 
The mar cs models also ma ke use of a slightly different equation 
of state ( Gustaf ssonl 1973 . and subsequent updates; hereafter 
"Uppsala EO S "). C ompared with the equation of state from 
Mihala s et al. (1988), the Uppsala EOS predicts for the most 



part slightly lower values of the gas pressure at a given gas 
density and temperature. 

The effects of the particular choice of equation of state on 
the calculation of the opacity tables are altogether negligible: 
differences in terms of bin opacities computed with the two 
equations of state are at most 0.5 % at solar metallicity and less 
than 0.1 % at [Fe/H] = -3 (for log 7-5000 > -4)- Our tests also 
indicate that the differences between the 3D-1D LTE correc- 
tions for Fe and O abundances determined from Fe 1 and OH 



lines using the equation of state from Mih alas et al.l 11988) and 
the Uppsala EOS are $0.02 dex for lines with equivalent widths 
smaller than 80 mA. 

Figure Q] shows the temperature structures as a function 
of optical depth resulting from the convection simulations at 
T e s — 4750 K, compared with the ID stratification from the 
corresponding marcs models. Figure [2] shows instead the tem- 
perature stratifications as a function of gas density for both 
the 3D simulations and ID marcs model atmospheres. At so- 
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Fig. 1. Thermal structures of four snapshots of 3D hydrodynamical simulations of red giants at different metallicities. Gray area: 
temperature distribution as a function of optical depth at X = 5000 A in the 3D convection simulations. Darker areas indicate 
temperature values with higher probability. Thin dashed line: Mean temperature stratifications of the 3D simulations (averaged 
over surfaces of equal optical depth at X = 5000 A). Thick solid line: Temperature stratifications of the corresponding marcs 
models. 



lar metallicities and mild metal-deficiencies ([Fe/H] k, -1), 
the mean temperature-density stratifications at a given opti- 
cal depth in the upper atmospheric layers of the hydrodynam- 
ical simulations closely resemble the ID structures of the cor- 
responding marcs models where radiative equilibrium is en- 
forced. At lower metallicities ([Fe/H] $ -2), instead, the tem- 
perature stratification in the outer layers of the simulations 
tends to remain significantly lower than in ID model atmo- 
spheres. The temperature in the optically thin layers of the con- 
vective simulations is for the most part regulated by two com- 
peting mechanisms: radiative heating caused by reabsorption 
by spectral lines of photons released at deeper layers, and adi- 
abatic cooling following the expansi on of the ascending gas. 



Following the interpretation given bv lAsplund etall d 19991) for 



3D hydrodynamical models of metal-poor solar-type stars, with 
fewer and weaker lines available at low metallicities, adiabatic 
cooling becomes more dominant and the balance between cool- 
ing and heating occurs at lower surface temperatures than in 
radiative equilibrium conditions. At [Fe/H] = -3, the average 
temperature difference between 3D and ID models in the up- 



per atmospheric layers (log T5000 ~ _ 3) is substantial and can 
amount to 1000 K or more at a given optical depth. 

The prediction of photospheric temperatures significantly 
below radiative equilibrium values is a crucial result of the 
convection simulations. In order to test the accuracy of the 
opacity binning scheme in this respect, we computed the ra- 
diative heating rates in "1.5D" approximation for all columns 
in vertical slices of four red giant simulation snapshots with 
varying metallicities using the opacity binning approach and 
monochromatic radiative transfer with opacity distribution 
functions (ODFs). The results are compared in Fig. [3] At very 
low metallicities ([Fe/H] S -2) there is an excellent correla- 
tion between the rates computed with the two approaches, sug- 
gesting that the opacity binning approximation is accurately re- 
producing the radiative heating and therefore also the tempera- 
tures at the surface of 3D metal-poor models. At solar metallic - 
ity and mild metal-deficiencies ([Fe/H] k, -1), the correlation 
is somewhat poorer, and the results plotted in Fig. [3] show that 
the binning scheme generally leads to underestimate both the 
heating and the cooling rates compared with monochromatic 
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Fig. 2. Thermal structures of four snapshots of 3D hydrodynamical simulations of red giants at different metallicities. Thin solid 
line: extreme temperatures at a given density in the 3D convection simulation. Thin dashed line: Mean temperature-density 
stratifications of the 3D simulations (averaged over surfaces of equal optical depth at A — 5000 A). Thick solid line: Temperature- 
density stratification for the corresponding marcs models. 



radiative transfer using ODFs. Such weaker response to ther- 
mal fluctuations with the opacity binning approach corresponds 
to a looser coupling between matter and radiation which in turn 
increases the relative importance of convective cooling. This 
possibly implies that, at these metallicities, the opacity binning 
scheme is likely to slightly underestimate the temperatures in 
the upper photospheric layers. 

Differences between the mean 3D and ID atmospheric ther- 
mal structures, as well as the presence of temperature and den- 
sity inhomogeneities in the 3D hydrodynamical simulations, 
can have dramatic effects on the predicted strengths of spec- 
tral lines. The cooler surface layers encountered in the con- 
vection simulations of metal-poor stars are expected to have a 
significant impact on temperature sensitive features, This is the 
case of, in particular, molecular lines and weak low-excitation 
lines of neutral metals, the line formation regions of which 
are shifted outwards in metal-poor hydrodynamical models be- 
cause of the lower photospheric temperatures. Also, as a con- 
sequence of the cooler mean temperature stratifications, in the 
upper photospheric layers of metal-poor stars, the gas and elec- 
tron pressure resulting from the convection simulations tend to 



rem ain lower than in the correspo nding ID hydrostatic models 
(see lAsplund & Garcia Pereal200ll) . The lower gas and elec- 
tron pressures in the metal-poor 3D hydrodynamical simula- 
tions are therefore expected to influence the formation of grav- 
ity sensitive lines. 

Figure 2] shows the spatially resolved outgoing intensity in 
the continuum bin for four snapshots of the simulations of the 
r e ff-4750 K series. The granulation pattern is clearly visible 
and qualitatively resembles the one observed on the Sun and 
in other simulations of late-type stars. The characteristic size 
of the granules varies depending on the metallicity of the sim- 
ulations. The size of granules may be shown t o scale approx- 
imately with the pressure scale height Hp (e.g. ISchwarzsc hild 
19751) . and with the ratio uh[u t of horizontal to vertical flow 
velocities dStein & Nordlund[l998t) . At low [Fe/H], the con- 
tinuous opacity is lower and thus the gas density is higher at 
a given optical depth. This means that smaller vertical veloc- 
ities are sufficient to sustain the convective flux. With smaller 
vertical velocities, horizontal velocities are also expected to be 
smaller. In practice, we observe that, in the proximity of the 
optical surface and below it (logTsooo ^ -1), the decrease in 
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Fig. 3. Comparison between the radiative heating rates computed at all grid-points in four vertical slices of four red giant sim- 
ulation snapshots with different metallicities using the opacity binning scheme ("BINS") and monochromatic radiative transfer 
with opacity distribution functions ("ODFs"). The radiative heating rates per unit mass (£2 m d) are normalized with respect to the 
specific heat (per unit mass) at constant pressure (Cp). The small boxes contain magnified views of the plots in the regions of low 
radiative heating rates. 



u z is overcompensated by a more pronounced decrease in uh, 
so that the product HpUn/u- becomes in fact systematically 
larger the higher the metallicity of the simulation, suggesting 
that granules must also be bigger. We also find that the size 
of the granules increases with the effective temperature of the 
simulations. We intend to return to a more detailed description 
of the physical properties of the convection simulations in a 
future paper. 



2.2. Spectral line formation 

We use the red giant convection simulations as time- 
dependent 3D hydrodynamical model atmospheres to per 
form detailed spectral line formation calculations. We fol 
low here the same procedure adopted in other recent in- 
vestigations of the effects of surface convect ion and gran- 
ulation on so lar and stellar spectroscopy ([A splund et al 



19991 l2000allbl I200I lAsplundl2000t lAsphind & Garcia Perez 
200U lAllende Prieto et al] l200ll l2002bt iNissen et al] |200 



lect representative sequences, typically ~100 to ~250 hours 
long (stellar time), of about 30 snapshots separated at regu- 
lar intervals in time. Prior to the line formation calculations, 
we decrease the horizontal resolution of the simulations from 
100x100 down to 50x50 to ease the computational burden. We 
also interpolate the simulations to a finer depth-scale, increas- 
ing the vertical resolution of the layers with log 7"5ooo^2.5 to 
improve the numerical accuracy. 



We compute flux profiles for lines of a number of ions 
and molecules under the assumption of LTE. As the main 
purpose of our study is to isolate and investigate the impact 
of 3D models on LTE spectral line formation, we have not 
performed calculations on a large selection of lines ordinar- 
ily used in abundance analyses. Instead we consider a sam- 
ple of "fictitious" atomic (Nai, Mgi, Cai, Fei, and Fen) and 
molecular (CH, NH, and OH) lines at selected wavelengths, 
with varying lower-level excitation potentials and line strength s 
dSteffen & Holwegerll2002t lAsplundll2005l: ICollet et al]|2006h . 



Collet et al.ll2006l) . From the full red giant simulations, we se- 



Concerning molecular lines, we restrict our investigation to 
a set of low-excitation (0 to 0.5 eV) features representa- 
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Fig. 4. Spatially resolved emergent intensity in the continuum bin for four snapshots of 3D hydrodynamical simulations of red 
giants at different metallicities; the characteristic surface granulation pattern is shown. In order to facilitate the comparison of 
the four cases, the patterns have been partially reproduced periodically so that the physical dimensions of the four images are the 
same as for the [Fe/H] = 0.0 simulation. 



tive of molecular bands frequently used in abundance anal- 
yses of giants (e.g. Ichristlieb et aill2004i iBessell et all l2004; 
Cavrel et all 12004: ISpite et alj|2005l) : CH lines at 4360 A be- 
longing to the A-X electronic transition band, the OH A-X 
system at 3150 A, and the NH A-X lines at 3360 A. Fictitious 
lines provide a bench-mark for a systematic comparison of ID 
and 3D LTE spectral line formation for various elements and 
molecules at different metallicities. They allow us to analyse 
the behaviour of spectral lines in ID and 3D models uniquely 
as a function of lower-level excitation potential, wavelength 
and line strength, separating it from other complications such 
as blends or wavelength dependence of continuous opacities. 
Finally, in addition to the fictitious spectral lines, we also con- 
sider some real transitions of particular interest for stellar spec- 
troscopy: the Li i line at 6707.8 A and the forbidden [Oi] lines 
at 6300.3 A and 6363.7 A. 



When computing ionization and molecular equilibria and 
continuous opacities for the line formation calculations, we as- 
sume the same chemical compositions as used for the construc- 
tion of the model atmospheres. As a rule, only the abundances 
of the trace elements are varied when calculating line opaci- 
ties for the 3D cases. The number densities of line absorbers 
are then calculated from Saha ionization and Boltzmann exci- 
tation balances and instantaneous molecular equilibrium at the 
local temperat ure. Partition function data for atoms and ions 
are taken from Irwinl ( 198lh and for molecu l es, tog ether with 
equilibrium constants, from lSauval & Tatuml d 19840 . 

The source function for lines and continuum is approxi- 
mated with the Planck function (S v — B v ) and scattering is 
treated as true absorption. We discuss the validity of this ap- 
proximation for the present analysis in Sect. 14.11 The radiative 
transfer equation is solved numerically along nine directions 
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(two //-angles, four </>-angles plus the vertical), after which we 
perform a disk integration and a time average over all selected 
snapshots. Various test calculations using a larger number of 
rays (four //-angles, eight ^-angles plus the vertical) indicate 
that, for the scope of the present analysis, our procedure is ac- 
curate enough in reproducing the spatially and temporally av- 
eraged line profiles; the differences between elemental abun- 
dances determined with our procedure and the test cases are 
typically less than 0.01 dex. In terms of derived abundances, 
temporal averaging over the selected number of snapshots is 
sufficient to obtain results with the same degree of accuracy, as 
verified by test calculations using different numbers of simula- 
tion snapshots. 

To estimate the impact of 3D hydrodynamical models on 
stellar spectroscopy, we perform a differential abundance anal- 
ysis using a simple curve-of-growth method. First, we calcu- 
late the LTE equivalent widths of the lines of our sample for a 
given chemical composition using ID marcs models. We then 
evaluate the 3D-1D LTE abundance corrections by varying 
the abundances in the 3D line formation calculations until the 
equivalent widths reproduce the ones computed in ID. Spectral 
line profiles are calculated for typically 60 to 100 wavelength 
points depending on the strength of the lines. Test calculations 
performed increasing the spectral resolution of the line profiles 
ensure that the accuracy of the differential 3D-1D analysis is 
consistently better than 0.01 dex with the adopted number of 
wavelength points. 

Particular care is exerted when dealing with OH and CH 
lines. The LTE number densities of the two hydrides in the up- 
per photospheric layers show a highly non-linear dependence 
not only on temperature but also on the relative abundances of 
carbon and oxygen. Because of the relatively large dissociation 
energy of the CO molecule (~11 eV), the formation of carbon 
monoxide can significantly reduce the densities of both oxy- 
gen and carbon available for OH and CH. Therefore, in order 
to properly evaluate 3D-1D effects for CH and OH lines, we 
are forced to take into account the simultaneous variation of 
carbon and oxygen abundances in the 3D line formation cal- 
culations. Here, we determine the 3D-1D corrections to C and 
O abundances self-consistently, by means of an iterative proce- 
dure. The analysis of NH lines, on the contrary, is not affected 
by such complications: changes in C and O abundances have, 
in fact, altogether negligible effects on the strength of NH lines. 
This simplifies the task of determining of 3D-1D corrections 
from these lines as variations of the N abundance can be stud- 
ied independently from the exact value of the C and O abun- 
dances. Vice versa, varying the N abundance in the spectral 
line formation calculations causes no appreciable changes on 
the strengths of CH and OH lines. 

The same numerical code and input physics are adopted for 
the line formation calculations with both ID and 3D model at- 
mospheres. In the ID cases, we compute spectral line profiles 
for two different values of the micro-turbulence, £ = 1.5 km s _1 
and £ = 2.0 km s _1 , corresponding to the typical range of 
values adopted in ordinary ID abundance analyses of red gi- 
ants with similar stellar parameters as the ones adopted for 
our suites of models. We emphasize, however, that none of the 
fudge parameters needed in the classical ID analyses (i.e. mix- 



ing length parameters, micro- and macro-turbulence) enters the 
3D spectral line formation calculations: only the velocity fields 
inherent to the hydrodynamical simulations are here used to 
reproduce non-thermal line broadening and asymmetries asso- 
ciated with convective Doppler shifts. 

For the fictitious lines we implement collisional broaden- 
ing by hyd r ogen atoms according to the classical description 
by lUnsoldl dl955 ) with an enhancement factor of 2.0 for the 
broadening constant. For the real features analysed here (Li i 
and [O i] lin es) we apply instead the quantum mechanical cal- 
culations by lBarklem et all feOOOh . 

As we limit our investigation to a differential analysis, we 
do not attempt to determine absolute 3D LTE abundances, but 
we rather address the question whether systematic effects are 
present in classical stellar spectroscopy based on ID models. 
The main advantage of performing a differential study is that 
the sensitivity of the results to uncertainties in background con- 
tinuous opacities, input physics, absolute line parameters (e.g. 
lower-level excitation potentials and log gf values), collisional 
broadening, or blends can be minimized. 

3. Results 

3.1. Fei and Fen lines 

Figure [5] shows the 3D-1D LTE corrections to the Fe abun- 
dance as derived, for the two series of model atmospheres, from 
Fe i lines at A = 5000 A with varying equivalent widths Wa- At 
a given wavelength, the magnitudes of the corrections and their 
trends as a function of line strength depend on several param- 
eters such as lower-level line excitation potential, metallicity, 
and effective temperature of the models. 

At a given Fe abundance, weak low-excitation Fe i lines ap- 
pear stronger in the framework of 3D models than they do in 
ID, resulting in negative 3D-1D LTE abundance corrections. 
The effects are more pronounced at lower metallicities where 
the differences between the 3D and ID Fe abundances can 

reach 1.0 dex. The 3D-1D abundance corrections become 

progressively smaller for higher-excitation lines, vanishing for 
line excitation potentials of approximately 5 eV. A similar re- 
sult characterizing granulation abundance co rrections for the 
Sun is reported by lSteffen & Holwe ger ( 2002). In t heir spectro 



scopic analysis of the red giant Pollux (ft Gem), Rula nd et al 



(1980) find that abundances derived for low-excitation lines 
are typically lower than for high-excitation lines and interpret 
the result as an effect of departures from LTE. Applying our 
3D- ID corrections, the trend of the derived abundances as a 
function of excitation potential would become even steeper, 
suggesting that the non-LTE effects might be even larger for 
the 3D case. Possible departures from LTE of Fe are discussed 
further in Sect. 1421 

The overall behaviour of the 3D- ID abundance corrections 
as a function line strength varies significantly depending on the 
metallicity of the models. At solar and moderately low metal- 
licities ([Fe/H] -1), the differences between 3D and ID 
abundances become progressively less negative with increasing 
line strength and may eventually revert sign, growing positive 
for sufficiently strong lines and relatively high excitation poten- 
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tials. It is also worthwhile observing that the abundance correc- 
tions for strong lines can reach a maximum and then turn-over 
as the line excitation potential increases. At lower metallicities, 
the 3D-1D abundance corrections generally show a shallower 
trend with equivalent width. In particular, at [Fe/H] = -3, 
the corrections can also grow more negative for stronger lines. 
The values of the 3D-1D abundance corrections for relatively 
strong lines (Wa ~ 50 mA) also depend on the choice of micro- 
turbulence, which controls the saturation level of the lines in 
the ID calculations. At a given abundance in ID the larger 
the micro-turbulence, the stronger the lines and in turn the less 
negative (or more positive) the differences between 3D and ID 
abundances become. 

Figure [6] shows analogous 3D-1D Fe abundance correc- 
tions computed for Fe i lines at 3500 A. The results are fairly 
similar to the ones for Fe i lines at 5000 A, the trends of the cor- 
rections following qualitatively the same behaviour as a func- 
tion of line strength. The exact values of the corrections can 
vary slightly for the two wavelengths due to the different steep- 
ness of the source function with depth and different sensitivity 
of the continuous opacity at 3500 A and 5000 A to temperature 
and electron density. 

In general, identifying the leading reason for a particular 
behaviour of the 3D-1D corrections as a function of the line or 
the model parameters is not a simple task, due to the complex- 
ity and non-linearities inherent to line formation, in particular 
in 3D hydrodynamical models. In practice, spectral lines form 
over a relatively large range of optical depths; in the frame- 
work of 3D models line formation is sensitive not only to the 
mean thermal stratification but also to inhomogeneities in the 
temperature, density, and velocity fields. The large abundance 
corrections for weak lines in the very metal-poor cases can be 
qualitatively understood by comparing the 3D and ID model 
structures. The mean temperature stratifications at the surface 
of 3D hydrodynamical simulations of very metal-poor giants 
are significantly cooler than in ID models computed for the 
same stellar parameters (Sect. 12. 11 1. This implies that in LTE 
the fraction of neutral to ionized iron is overall enhanced in 
the upper layers of 3D metal-poor models compared with the 
ID case. The cooler temperature structure of 3D models also 
contributes to reducing the electron pressure and, consequently, 
lowering the density of H ions, thus decreasing the opacity in 
the continuum-forming layers. Therefore, at low metallicities, 
the combined effect of increased line opacity and decreased 
continuous opacities tends to make weak low-excitation Fei 
lines stronger in 3D than in ID; a lower Fe abundance is then 
required in metal -poor 3D models to reproduce the ID equiva- 
lent widths. 

Weak high-excitation lines form in deeper photospheric 
layers and are less sensitive to the temperature structure at the 
surface of the models. At all metallicities, the resultant 3D-1D 
abundance corrections for high-excitation lines are smaller than 
for the low-excitation features. Stronger lines form higher up in 
the atmosphere compared with weak lines when the other line 
properties remain the same. While at very low metallicities Fe n 
is clearly the dominant Fe ionization stage nearly everywhere 
in the ID models, in 3D the Fei fraction is substantial at all 
depths contributing to the emergent line profiles. Therefore in 



very metal-poor models the 3D-1D Fe abundance corrections 
typically remain negative for strong Fei lines. As a test, we 
compute curves-of-growth of Fe i lines using the horizontally 
averaged 3D structure to quantitatively study the dependence of 
line strengths on the lower temperatures of the mean stratifica- 
tion. The behaviour of the 3D-1D corrections for Fe i lines can 
be qualitatively reproduced by comparing the curves-of-growth 
for the ID marcs and mean 3D stratifications. The test confirms 
that the lower temperatures encountered in the 3D hydrody- 
namical metal-poor model atmospheres are the main factor de- 
termining the large and negative 3D-1D corrections. However, 
we emphasize that, from a quantitative point of view, line for- 
mation calculations relying on the mean 3D stratification can- 
not accurately reproduce the results of the calculations based 
on the actual 3D model atmospheres. In fact, abundances de- 
rived from weak low excitation Fe i lines using the full 3D and 
the mean 3D stratifications can differ by up to about 0.3 dex. 
This is an indication that the temperature and velocity inhomo- 
geneities of the 3D structure should not be neglected. 

At metallicities closer to solar, the trends of 3D-1D Fe 
abundance corrections with line strength cannot be ascribed 
simply to differences between ID and mean 3D stratifications, 
which are very similar in the line forming layers. Curves-of- 
growth computed for Fe i lines using either the ID marcs mod- 
els or the mean 3D stratifications (and the same value of micro- 
turbulence) closely resemble each other. The differences in Fe 
abundance determinations using the ID and mean 3D strati- 
fications amount in fact to less than 0.04 dex for Fei lines 
within the range of line strengths considered here. Hence, the 
actual trends of the 3D- ID corrections at metallicities near 
solar must be attributed mainly to 3D temperature inhomo- 
geneities and possibly velocity fields. We also note that the 
trend of the 3D- ID abundance corrections, increasing with 
equivalent width for [Fe/H]S - 1, can be flattened by appropri- 
ately choosing a smaller value for the micro-turbulence than the 
ones considered here (£ < 1 .5 km s ). However, typical micro- 
turbulences derived from spectroscopic ID analyses of red gi- 
ants are found indeed to be in the range between 1 .5 km s _1 and 
2.0 km s -1 . At metallicities [Fe/H]£ - 1, the present 3D-1D 
Fe abundance corrections, would therefore introduce a trend in 
the derived abundances as a function of equivalent width. This 
could be an indication that convective motions of the gas are 
not fully resolved in the present hydrodynamical simulations 
and, consequently, that non-thermal Doppler broadening is also 
underestimated. Further investigation is however necessary in 
order to identify other possible causes for these trends. 

We observe that the 3D- ID Fe abundance corrections for 
lines at a given equivalent width are typically shifted toward 
more negative values for the higher effective temperature mod- 
els. The equivalent width of a (weak) spectral line is roughly 
proportional to the ratio { y /k v of fine to c ontinuous opacity at 
line formation depth (see also Gravll 19921 p. 277). In the case 
of (weak) low-excitation Fe i lines, for a given Fe abundance, 
the ratio (£ v /*v)3d / (^v/*v)id turns out to be larger the higher the 
effective temperature of the models at all metallicities, indicat- 
ing that 3D-1D Fe abundance corrections are also expected to 
be more pronounced (more negative). 
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The 3D- ID Fe abundance corrections for Fen lines at 
5000 A and 3500 A are presented in Fig.|7]and[8] The overall 
trends of the corrections for these lines are similar at all metal- 
licities; differences between abundances derived with 3D and 
ID models are typically positive and generally increase with in- 
creasing line strength and excitation potential. Low-excitation 
lines can, however, show negative corrections at metallicities 
below [Fe/H] -2. Similarly as for Fe i lines also, the correc- 
tions Fe ii lines at a given equivalent width can reach a maxi- 
mum and revert their trend for sufficiently high excitation po- 
tentials. This behaviour is also present in 3D- ID abundance 
corre ctions for lines o f neutral and ionized metals in solar-type 
stars dAsplundll2005l) . although evident already for somewhat 
weaker lines. 



Table 2. Differential 3D and ID Li abundances derived from 
the Li i 6707.8 A line. The computed equivalent widths Wa of 
the Li i line are also repo rted. A log gf value of 0.174 is as- 
sumed |Smiftetal] {l99a) in the line formation calculations. 



(T eS ) 


[Fe/H] 


loge(Li) 1D 


log e(Li) 3D 


W A 


[K] 








[mA] 


4697 


+0.0 


1.0 


0.90 


37.7 


4717 


-1.0 


1.0 


0.77 


34.7 


4732 


-2.0 


1.0 


0.62 


30.1 


4858 


-3.0 


1.0 


0.55 


22.0 


4983 


+0.0 


1.0 


0.78 


17.5 


5131 


-1.0 


1.0 


0.66 


11.7 


5035 


-2.0 


1.0 


0.56 


14.7 


5128 


-3.0 


1.0 


0.50 


12.7 



3.2. Nai, Mgi, and Cai 

Figures [9j[10] and [TTI show the 3D-1D LTE corrections to the 
sodium, magnesium, and calcium abundances as derived from 
Na i, Mg i, and Ca i lines at 5000 A. The behaviour of the abun- 
dance corrections as a function of line strength for these lines 
is similar to the one predicted for Fe i lines. In LTE, the over- 
all trends of the corrections are governed by the differences 
between the ID and mean 3D temperature stratifications and 
by the presence of temperature inhomogeneities and velocity 
fields in the 3D models. The actual magnitude of the 3D-1D 
effects, however, largely depends on the ionization potential of 
the species under consideration. The 3D- ID abundance cor- 
rections for Na i and Ca i, for instance, are considerably smaller 
than the ones derived from Fe i lines. In fact, the relatively low 
ionization potentials of Na i and Ca i imply that the two metals 
are more easily ionized compared with Fe i, both in ID and 3D 
models. This significantly reduces the 3D-1D effects on line 
strengths due to temperature inhomogeneities or differences in 
the thermal stratifications. Magnesium on the contrary has a 
fairly high first ionization potential, comparable to that of iron: 
hence, the computed 3D-1D abundance corrections for Mgi 
lines are typically also as large as the ones derived from Fe i 
lines. 



3.3. Li i 

Lithium abundances in stars are of great importance in as- 
trophysics: they can be used as diagnostics of stellar evolu- 
tion, primordial nucleosynthesis as well as cosmic and Galactic 
chemical evolution. Nearly all stellar Li abundance determi- 
nations in late-type stars are based on the analysis of the Lii 
resonance line at 6707.8 A. In order to accurately determine 
Li abundances it is therefore necessary to properly model the 
formation of this line, taking into account possible effects due 
to stellar granulation as well as departures from LTE in gen- 
eral. Recent investigations of 3D/1.5D non-LTE Lii spectral 
line f ormation in the Sun (|Kiselmanl 1 19971 119981: llJitenbroekl 
1998 | ) and metal-po or solar-type stars ( Asplund et al .1 120031: 
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20031) indicate that non-LTE effects on the de- 
rived abundances are significant and comparable to the ones 
due to granulation. In the present study, however, we only 



[Fe/H] 1D 

Fig. 12. 3D-1D LTE corrections to Li abundances derived 
from the Li i line at A = 6707.8 A as a function of metallic - 
ity [Fe/H] of the models. 



explore the effects of granulation on the formation of the 
Li i 6707.8 A line under the approximation of LTE; a detailed 
3D non-LTE analysis of Li i line formation in giant stars is de- 
ferred to a future paper. In Table [2] (Fig. IT2t . we present the 
3D-1D LTE corrections to Li abundances derived from the 
Lii 6707.8 A line, assuming a value of loge(Li) 1D = 1.0 at all 
metallicities in the ID calculations. The overall 3D-1D LTE 
Li abundance corrections as a function of metallicity behave 
rather similarly to the ones derived for weak (iS40 mA) low- 
excitation lines of other neutral elements. Since, the ioniza- 
tion potential of Li i is just slightly larger than the one of Na i, 
the 3D-1D LTE corrections for these two species are typically 
very similar. 

Once again, the differences between the 3D and ID abun- 
dances are relatively large and negative ($-0.5 dex) at very low 
metallicities, because of the markedly cooler temperature struc- 
tures of 3D model atmospheres compared with the ID cases, 
while they become gradually smaller as metallicity increases. 
Also, as for neutral lines in general, 3D-1D LTE Li abundance 
corrections are more pronounced for the series of model atmo- 
spheres with higher T e s. 
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Fig. 13. 3D-1D LTE corrections to O abundances derived from the forbidden [Oi] lines at A — 6300.3 A and A = 6363.7 A as a 
function of metallicity [Fe/H] of the models. 



Table 3. Comparison of the 3D and ID O abundances 
derived from the [Oi] forbidden lines at 6300.3 A and 
6363.7 A. The log g f va lues of the two lines are adopted from 
Storey & Zeippenl (120001) . The 3D O abundances are the ones 



that reproduce the equivalent widths computed with ID models 
for the tabulated ID abundances. 
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loge(0) 1D 


loge(0) 3D 


loge(0) 3D 
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6300.3 A 


6363.7 A 
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+0.0 


8.89 
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4858 


-3.0 


5.89 


5.77" 
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8.36 
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4732 
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7.33 


7.33 
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6.39 
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6.23 a 


5131 


-1.0 
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The equivalent width of the line is less than 1 mA. 



as the ones considered hereQ The existence of oxygen over- 
abundances with respe ct to iron in metal- poor stars has been 
pointed out already by IConti et al.1 (119671) . The exact amount 
of the over-abundances, however, and their overall trend with 



metal licity is still today a matter of debate dGarcfa Perez et al 



2006). In our analysis, we therefore consider, for simplicity, 



two different values of oxygen enhancement with respect to the 
scaled solar standard compositions for the ID spectral line for- 
mation calculations: [O/Fe] = +0.0 and [O/Fe] = +0.5 (for 
[Fe/H]< - 1). 

In Table |3] we present the differential 3D and ID O abun- 
dances as derived from the two [Oi] lines. The trends of the 
3D-1D abundance corrections as a function of metallicity or 
effective temperature of the models (Fig. [L31 are qualitatively 
the same as in the case of the other neutral species discussed 
above. Nonetheless, the 3D-1D corrections for [Oi] forbidden 
lines are in general significantly smaller. Because of its high 
ionization potential (13.6 eV), Oi is, in fact, the dominant ion- 
ization stage for oxygen in the line formation layers of both 
3D and ID models and shows rather little sensitivity to the dif- 
ferences between the 3D hydrodynamical and ID hydrostatic 
temperature stratifications at the surface (see also Sect. 13.21 1. 
The small differences between the two [O i] lines in terms of 
abundance corrections are due to the different formation depths 
of the two features. 



3.5. CH, NH, and OH lines 



3.4. [0\] lines 

The [Oi] forbidden lines at 6300.3 A and 6363.7 A are among 
the most widely used indicators of O abundance in c ool stars 
and halo giants (r e ff ~ 5000 K) in particular (e .g. Barbuy 



1988[lAllende Prieto et al.l200lHNissen et alkOOalSpite et al 
2005t iGarcfa Perez et al.l 120061) The main advantage concern- 
ing the use of [ O i] lines is that th ey are expected to form under 
LTE conditions ( Kiselmanl2001 ): the drawback is, on the other 
hand, that these lines tend to become, in general, very weak 
already at [Fe/H] ^ -2 in stars with similar surface gravities 



Molecule formation shows an extreme and highly non-linear 
sensitivity to temperature in the upper layers of late-type stel- 
lar photospheres. Because of this strong temperature depen- 
dence, molecular lines are very susceptible to systematic er- 
rors in the thermal structure of the model atmospheres. In 3D 
model atmospheres, the vertical temperature gradients are typ- 



1 At [Fe/H] = -2 and with [O/Fe] = 0, the equivalent widths of the 
two [Oi] lines at 6300.3 A and 6363.7 A computed using the marcs 
model atmosphere at T cS = 4732 K are 3.7 mA and 1.2 mA, respec- 
tively. 
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Fig. 14. 3D-1D LTE corrections to O and C abundances derived from fictitious OH lines at A = 3150 A and CH lines at 
A = 4360 A as a function of metallicity [Fe/H] of the models. 



ically steeper in upflows than in downdrafts. The temperature 
contrast reverses in the upper photosphere some distance above 
the optical surface, so that the gas above the granules becomes 
coole r than average in the overshoot region (iStein & Nordlund 
19981) . The temperature contrast in these layers is more pro- 
nounced for the lower metallicity models where the coupling 
between the radiation field and the gas is weaker. Temperatures 
in the high photospheric layers are generally low; however, 
compression due to converging gas flows or shocks at the 
boundaries between granules heat the fluid in the intergranu- 
lar lanes, raising the temperature above the radiative equilib- 
rium value. Consequently, the strengths of molecular lines are 
expected to vary across the surface of the model atmospheres, 
being stronger in the upflowing regions. Also, the formation of 
the spatially and temporally resolved molecular line profiles is 
biased towards granules because of their large area coverage, 
steeper temperature gradients, and higher continuum intensi- 
ties. 

The significantly lower temperatures of the upper photo- 
spheric regions of 3D metal-poor model atmospheres com- 
pared with their ID counterparts favour higher concentrations 
of CH, NH, and OH molecules for a given chemical compo- 
sition, therefore resulting in stronger molecular lines and neg- 
ative 3D-1D C, N, and O abundance corrections. Differences 



between the 3D and ID predicted molecular line strengths are 
less pronounced at higher metallicities, in accordance with the 
similarities between the ID and mean 3D temperature stratifi- 
cations, and are ascribable primarily to the temperature inho- 
mogeneities in the 3D structures. 

In Tables |4]and|5]we present the 3D-1D LTE C, N, and O 
abundance corrections derived from weak fictitious CH, NH, 
and OH lines (with equivalent widths W.1^30 mA), as a func- 
tion of metallicity (see also Fig.[T4land[T5ll. At very low metal- 
licities ([Fe/H] -3), the 3D- ID corrections are very large 
and negative and range between -0.5 and - 1 .2 dex for the cases 
here considered. As expected, the differences between 3D and 
ID abundances become gradually less pronounced at higher 
metallicities, reducing to values of approximately -0.1 dex at 
[Fe/H] = 0. The actual magnitude of the abundance corrections 
for molecular features depends on several factors such as the 
species under consideration, the line parameters (in particular 
the lower level excitation potential), and the effective temper- 
atures of the models. The 3D- ID effects for higher excitation 
lines at a given metallicity are generally smaller in accordance 
with the sensitivity of these lines to deeper layers where the 
temperature differences between 3D and ID models are less 
prominent. The overall larger corrections for hotter models are 
in agreement with the larger 3D-1D temperature differences in 



20 



R. Collet et al.: 3D Hydrodynamical Simulations of Surface Convection in Red Giants 



Table 4. Differential 3D and ID C and O abundances derived from the analysis of weak fictitious CH and OH lines (see 
Sect. |2.21 >. The 3D abundances are the ones that reproduce the equivalent widths of fictitious lines calculated with ID models for 
the tabulated ID C and O abundances. 
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very metal-poor stars and, at near solar metallicities, with the 
presence of temperature inhomogeneities in the 3D structure 
and the high non-linearity of molecule formation. 

The magnitudes of the corrections for CH, NH, and OH 
lines naturally depend on the details of the molecular equilib- 
rium with other species, and the various indicators in practice 
show different sensitivities not only to the photospheric tem- 
perature structure but also to the chemical composition. For in- 
stance, at very low metallicities, the corrections for NH and 
CH lines differ by up to 0.2 dex in spite of the similar dissoci- 
ation energies of the two molecules (3.42 eV and 3.47 eV, re- 
spectively). Line strengths of CH and OH molecules are also 
sensitive to both carbon and oxygen atomic abundances. In 
Table l4l (Fig. [T4li we compare the results of 3D- ID analyses 
on OH and CH lines performed assuming, for the ID calcula- 
tions, two different values of oxygen enhancement with respect 
to the scaled solar standard compositions: [O/Fe] = +0.0 and 
[O/Fe] = +0.5 (for [Fe/H]< - 1), with all other metals regu- 
larly scaled as usual by the same amount as iron. Changes in 
the oxygen abundance alter the number densities of of both OH 
and CH molecules by different amounts in 3D and ID. In par- 
ticular, increasing [O/Fe] from +0.0 to +0.5 in the ID calcu- 
lations results in overall larger 3D-1D abundance corrections 
for OH lines and smaller corrections for CH lines. In practice, 
however, the effects on the abundance corrections are appre- 
ciable only at very low metallicities. The trends of 3D-1D C 
and O abundance corrections from CH and OH molecular lines 
at low metallicity depend on the particular choice made here 
for the ID chemical compositions (with C/O < 1) and on the 
sensitivity of molecule formation to temperature in the coolest 
layers of 3D models. In the very metal-poor ID giant mod- 
els, most of the carbon and oxygen are in atomic form, due to 
the low surface gravity, the overall low C and O abundances, 
and the relatively high temperatures of the upper photospheric 
layers (~4000 K or more). The fraction of carbon and oxy- 
gen atoms locked in CO molecules is negligible and increases 



only slightly when [O/Fe] is changed from to +0.5. As a 
result, the ID photospheric number density of CH molecules 
also remains approximately the same in both the [O/Fe] = 
and [O/Fe] = +0.5 cases. On the contrary, in very metal-poor 
3D models, temperatures can be significantly lower in the up- 
per photosphere and, below ~3700 K, a considerable fraction 
of carbon ends up locked in CO. As [O/Fe] increases from to 
+0.5, the number density of CO in these cool layers increases 
further, reducing the amount of carbon available for the for- 
mation of CH molecules. Thus, while the number density of 
CH molecules is still larger than in ID, the 3D- ID C abun- 
dance corrections become smaller (less negative) if [O/Fe] in- 
creases. In the case of oxygen, in the very metal-poor ID mod- 
els where the fraction of CO is negligible, the number density 
of OH molecules in the photosphere simply scales proportion- 
ally to the O abundance, i.e. it increases by a factor of 10 5 ^3 
when [O/Fe] is raised from to +0.5. In the coolest layers 
of 3D models, where most of the carbon is locked in CO, the 
amount of oxygen available for OH molecule formation is also 
reduced. Assuming log(C/0)Q ^ -0.3 and that most carbon 
is in form of CO molecules, the amount of oxygen available to 
form OH molecules or remain in atomic form increases by a 
factor as high as (10 8 - 1)/(10 3 - 1) ^5 > 3 as [O/Fe] is 
changed from to +0.5; that also explains why the 3D- ID O 
abundance corrections derived from [0 1] and OH lines become 
more pronounced as [O/Fe] increases. 

4. Restrictions of line formation calculations 

4.1. Scattering 

The results of our analysis indicate that the differences between 
abundances derived from 3D and ID model atmospheres of red 
giant stars can be considerable, especially at very low metal- 
licities. The presence of large 3D- ID corrections can have 
strong implications for the interpretation of stellar abundances 
in terms of Galactic chemical evolution. The current 3D anal- 
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Table 5. Differential 3D and ID N abundances derived from 
the analysis of weak fictitious NH lines. The 3D abundances are 
the ones that reproduce the equivalent widths of fictitious lines 
calculated with ID models for the tabulated ID N abundances. 
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Fig. 15. 3D- ID LTE corrections to N abundances derived from 
fictitious NH lines at A = 3360 A as a function of metallicity 
[Fe/H] of the models. 



ysis, however, still relies on a number of approximations and 
assumptions, which can possibly lead to systematic errors in 
the derived abundances and must therefore be discussed. 

As mentioned in Sect. |2.1| and |2T2l we treat scattering as 
true absorption when solving the radiative transfer equation, 
both for the convection simulations and the line formation cal- 
culations with 3D and ID model atmospheres. This assumption 
can in general lead to systematic errors in the predicted tem- 
perature stratification of the upper layers of 3D hydrodynam- 
ical model atmospheres where the contribution of scattering 
to extinction is significant, and, ultimately, affect the strengths 
of synthetic spectral lines. Rayleigh scattering of Hi can con- 
tribute significantly to the total continuous extinction in the 
UV and blue part of the spectrum. At these wavelengths, im- 
plementing scattering as true absorption in the detailed line 
formation calculations underestimates, in general, the emer- 
gent flux in the con tinuum, resulting in too weak spectral lines 



(ICayrel et alJ 120041) . The effect is expected to be more pro- 



nounced in very metal-poor stars, where line-blocking in the 
continuum is weak, and, in particular, in metal-poor 3D model 
atmospheres, where the lower surface temperatures result in 
lower electron number densities and higher densities of scat- 



terers (H i particles) relative to absorbers (namely H particles 
in the case of continuous absorption). 

In order to quantify the effect of our approximated treat- 
ment of scattering on the photospheric stratifications, we have 
computed a series of ID marcs model atmospheres of giant 
stars at various metallicities, including scattering as true ab- 
sorption in the solution of the radiative transfer equation. The 
results of our test indicate that treating scattering as true ab- 
sorption leads to hotter temperature stratifications compared 
with mod els in which scattering i s properly taken into account 
(see also Gustafsson et al. 1975b . Temperature differences at 
optical depth log T5000 = - 3 can reach 300 K for model stel- 
lar atmospheres of very metal-poor red giants at [Fe/H] = -3; 
for [Fe/H] = — 1 model atmospheres, the temperature differ- 
ences at the same optical depth are smaller and amount to 
about 100 K. On the sole basis of these results, one could 
be tempted to conclude that the temperatures at the surface 
of the metal-poor 3D hydrodynamical models might be over- 
estimated. However, further testing is needed to estimate the 
magnitude of the effect on the 3D temperature stratification in 
metal-poor stars, and, more importantly, to also assess whether, 
in this case, it proceeds in the same direction. Either way, the 
treatment of scattering as true absorption is likely to signifi- 
cantly affect the balance between radiative heating and adia- 
batic cooling in the upper photospheric layers. 

We have also performed some test line formation calcula- 
tions with the average 3D stratifications of our very metal-poor 
model red giant stellar atmospheres and using the the spectrum 
synthesis code bsyn from the Uppsala package, which treats 
continuum scattering properly within the LTE approximation. 
We find that including scattering as true absorption in the cal- 
culations has relatively little influence on the predicted line 
strengths, at least for the type of giants considered in our study. 
Nonetheless, the effect of the above assumption might, in prac- 
tice, be larger in full 3D line formation calculations because of 
the presence of temperature and density inhomogeneities. The 
use of a differential 3D-1D abundance analysis ensures, how- 
ever, that the uncertainties in the treatment of scattering in the 
spectral line formation calculations are at least reduced. 

4.2. Departures from LTE 

It is important to emphasize that, in general, spectral lines of 
many of the species considered in our analysis are expected to 
suffer from departures from LTE. In particular, Fe 1 lines are 
most certainly seriously affected by non-LTE effects. The main 
non-LTE mechanism for Fe in late-type stellar atmospheres is 
over-ionization, which is driven by the radiation field J v in the 
UV being larger than the Planck function at the local temper- 
ature B V (T). This causes efficient photo-ionization from Fei, 
leading to underpopulation of all Fei levels compared with 
LTE and to weaker Fe 1 lines. This, in turn, implies that Fe 
abundances derived from Fe 1 lines are generally larger in non- 
LTE than with the assumption of LTE. At present, there is no 
consensus on how severe the non-LTE effects actual ly are on 



Fei l i nes in late-type star s (see, e.g., discussion in lAsplund 
20051) . ICollet et alJ (I2006T) have estimated the non-LTE effects 
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on Fe i lines for the extremely metal -poor giant HE 0107-5240 
by means of a ID analysis based both on a marcs model atmo- 
sphere and the mean atmospheric stratifications inferred from 
3D simulations. The non-LTE effects there are found to be con- 
siderable and opposite to the 3D-1D LTE corrections, confirm- 
ing that a combined treatment of both 3D and non-LTE effects 
is necessary for accurate Fe abundance determinations. A full 
3D non-LTE study of Fe i spectral line formation is certainly 
of high priority, but goes beyond the scope of the present work 
and is postponed to a future paper. 

Similarly, lines of other species considered in this study, 
such as Li i, Na i, Ca i, and Mg i might form out of LTE con- 
ditions. Molecular lines might as well suffer from departures 
from LTE; however, non-LTE effects on molecular line for- 
mation are, still today, largely u nexplored even in ID analyses 
dAsplund & Garcia Perezll200ll) . In addition, molecule forma- 
tion could also occur out of LTE: the rapid cooling occurring 
at the surface of the 3D simulations could imply that molecular 
equilibrium is not reached; also, photo-dissociations feeding on 
the radiation field coming from deeper layers might move the 
molecular equilibrium out of LTE. These aspects are undoubt- 
edly worthy further investigation. 



4.3. Stellar parameters 

The derivation of accurate abundances requires not only re- 
alistic modelling of the structure of stellar atmospheres and 
of the processes of spectral line formation, but also appro- 
priate determination of the fundamental stellar parameters. In 
the differential 3D-1D analysis presented here, we have as- 
sumed the stellar parameters to be the same for both 3D and 
ID model atmospheres. However, one could expect the deter- 
mination of the stellar parameters for a given star to depend 
on whether ID or 3D models are used (see also discussion in 
Asplund & Garcia PerezfcoOl). For instance, the results of our 
differential analysis of Fe lines suggest that differences in the 
derived metallicities cannot indeed be excluded. Also, differ- 
ences between the thermal structures of 3D and ID model at- 
mospheres can affect ionization equilibria and, therefore, the 
determination of spectroscopic surface gravities. Finally, tem- 
perature inhomogeneities in the continuum-forming layers of 
3D hydrodynamical model atmospheres can, in general, lead to 
different emergent flux distributions compared with the predic- 
tions of ID models and, consequently, affect effective temper- 
ature estimates. 



The above discussion on the possible sources of systematic 
errors in the 3D modelling of stellar atmospheres and line for- 
mation can create the impression that the results presented here 
are rather uncertain and that analyses based on ID model atmo- 
spheres are still preferable. It is important, however, to empha- 
size that many of the uncertainties in the present 3D analysis 
(e.g. departures from non-LTE and molecular equilibrium) ap- 
ply as well to the ID analysis, besides the systematic errors 
introduced, in the latter case, by the assumption of hydrostatic 
equilibrium and the rudimentary treatment of convective en- 
ergy transport. 



5. Conclusions 

In the present work, we have investigated the impact of 3D 
hydrodynamical model atmospheres of red giant stars on the 
formation of spectral lines of various ions and molecules and 
on the derivation of elemental abundances. The differences be- 
tween the mean 3D temperature stratifications and correspond- 
ing ID model atmospheres as well as the presence of tempera- 
ture and density inhomogeneities and correlated velocity fields 
in the 3D structure, can significantly affect the predicted line 
strengths and, hence, the values of elemental abundances de- 
rived from spectral lines. Differences between abundance de- 
terminations based on 3D and ID models are particularly large 
at low metallicities. The temperatures in the upper layers of 
3D hydrodynamical model atmospheres of metal-poor giant 
stars are in fact significantly lower than the ones predicted by 
classical ID models computed for the same stellar parameters. 
Because of the cooler temperature stratifications of 3D model 
atmospheres, weak lines of neutral species and molecules ap- 
pear considerably stronger in 3D than in ID, under the assump- 
tion of LTE. The 3D-1D corrections for these lines are large 
and negative. In particular, we find the 3D-1D LTE correc- 
tions to CNO abundances derived from CH, NH, and OH weak 
low-excitation to be typically in the range -0.5 dex to - 1 .0 dex 
for the very metal-poor giants at [Fe/H] ^ -3 considered here. 
We also derive large negative corrections to the Fe abundances 
(about -0.8 dex) at these metallicities, based on the analysis of 
weak low-excitation Fe i lines. 

In Sect. 2J we finally discuss possible systematic errors af- 
fecting the present 3D abundance analyses. A major source 
of uncertainty in this respect is the approximated treatment of 
scattering as true absorption, both in the 3D convection sim- 
ulations and in the line formation calculations. As discussed 
in Sect. 14.11 this approximation can lead to systematic errors 
in the estimated temperatures of the upper layers of 3D model 
atmospheres as well as in the predicted emitted fluxes in the 
UV, therefore the strengths of spectral lines. Also, the assump- 
tion of LTE is in general questionable for many of the ions and 
molecules considered in the analysis. In particular, in the case 
of iron, strong over-ionization feeding on the UV radiation field 
can lead to significantly weaker Fe i lines in metal-poor red gi- 
ants stars compared with the cases where LTE is enforced. The 
effects of over-ionization on the strengths of Fe I lines are op- 
posite to the ones of granulation. Departures from LTE should, 
therefore, accounted for in the 3D line formation calculations 
in order to accurately determine 3D-1D corrections to the Fe 
abundances. 
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